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ABSTRACT 

We compute to next-to-leading order in QCD the tensor unparticle contribution to 
the diphoton production at the LHC, wherein the unparticle sector is a consequence 
of (a) scale invariance but not full conformal invariance and (b) conformal invariance. 
We use the semi-analytical two cutoff phase space slicing method to handle the 0(a s ) 
corrections to the p p — > 7 7 X and show that our results are insensitive to the soft and 
collinear cutoffs. In order to avoid the contribution of the photons due to fragmentation, 
we employ the smooth cone isolation criterion. Significance of the QCD corrections to 
the diphoton events including unparticles is highlighted. 
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1 Introduction 



Diphoton signals play very important role in the phenomenological studies of the stan- 
dard model (SM) and also of beyond the standard model (BSM). These photons can 
constitute a potential background to those coming from Higgs decay. In the BSM, 
namely, supersymmetric [1], extra-dimensional [2] and unparticle (UP) models [3] the 
diphoton productions can constrain the model parameters. Hence it is important to 
have precise predictions [4-6] for this process in the theory in order to discover the Higgs 
boson as well as to unravel physics beyond the standard model. Recently, unparticle 
physics [7, 8] has gained a lot of attention in the context of collider signatures aiming to 
study new physics signals. In a recent paper [3], it was found that the diphoton pro- 
duction at the hadron colliders is an important process to look for effects coming from 
unparticle physics. This study was based on a leading order (LO) computation in QCD. 
Since the gluons contribute even at LO in the unparticle scenario, the higher order QCD 
effects are expected to play a significant role in the predictions. In this paper we report 
contributions to next-to-leading order (NLO) in QCD, coming from a model with tensor 
unparticles as well as from the SM. 

Parton level cross sections beyond the LO in perturbation theory encounter IR sin- 
gularities resulting from the soft gluons and massless collinear partons. In addition to 
the above singularities resulting from partons, the photons in the final states can also 
give collinear singularities when they become collinear to one of the final state par- 
tons. These collinear singularities, often called QED singularities, can be avoided by 
appropriately constructing observables that are insensitive to them. In addition, a par- 
ton in the final state can fragment non-perturbatively into a photon and a bunch of 
hadrons collinear to it, parameterized by photon fragmentation functions [9]. The QED 
singularities can be removed by either absorbing them into these fragmentation func- 
tions or by suitable definition of diphoton events that removes these contributions. The 
approach involving fragmentation functions is theoretically tedious and also brings in 
additional non-perturbative input that is poorly known to date. Hence we adopt an 
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alternate smooth cone isolation criterion proposed by Frixione [10] which ensures that 
the fragmentation contribution and the final state QED singularities are suppressed. 

2 Unparticle physics 

In a recent paper, Georgi [7, 8] considered a scenario wherein the SM is weakly coupled 
to a sector which is scale invariant below an IR scale Au- Consequently conventional 
particle interpretation at low energies in this new sector fails. These fields were termed 
as unparticles by Georgi. The unparticles could be interpreted as fractional number of 
massless particles [7] or a collection of particles with a particular mass distribution [11] 
as a result of deconstruction. The coupling of these unparticles to the SM would lead 
to direct phenomenological consequences due to the peculiar unparticle phase space [7] 
and propagator [8, 13], which are determined by scale invariance. Their coupling to the 
SM Higgs could lead to a constraint on the unparticle sector due to breaking of the 
scale invariance when the Higgs gets a vacuum expectation value [12]. As a consequence 
unparticle physics may be relevant at the high energy colliders which have been actively 
investigated [3,7,8,13-15]. In [16], the corrected propagator and lower bounds on the 
scaling dimensions of the vector and the tensor operators were obtained using unitarity 
constraints on scattering amplitudes. In this paper, we use these new corrected results 
to do our phenomenology. 

The unparticle operators could be of scalar, vector, tensor or fermionic type. In the 
present article we restrict ourselves to tensor unparticles [3,14] coupled to SM fields given 
by 

y^ ^' (1) 

where A t is the dimensionless coupling constant for the unparticle tensor operator Off 
which is traceless and symmetric and has a scaling dimension du- T^ v is the energy 
momentum tensor of the SM. Scale invariance restricts the scaling dimension of tensor 
unparticle operator to du > 3 [16]. Conformal invariance on the other hand leads to a 
constraint du > 4 on the second rank tensor operators. Scale and conformal symmetries 



can only guide us on fixing the tensor structures of the propagator leaving the overall 
normalisation undetermined. Unlike the conformal invariance, the scale invariance does 
not fix the relative coefficients of the tensors appearing in the tensor propagator [16]. 
We use the following propagator for our phenomenological study: 

Jd*xe^(0\TOi?(x)Of(0)\0) = ~iC T ^~ {d ^ 2) (-k 2 )^ 2 

x [du(d u - l){g m 9vi3 +fi++u) + ...] . (2) 

where we have chosen Ct =1. The terms given by ellipses do not contribute to the 
diphoton production. The terms in the ellipses depend on tensors, proportional to g^ v , k 11 
and k v . The exact tensorial form of course depends on the symmetry (scale or conformal). 
These terms do not contribute to physical processes thanks to the conservation and 
traceless nature of the SM energy momentum tensor. Hence the symmetry restriction 
enters only through the scaling dimension du (the overall undetermined constant could 
be different for the scale and conformal invariant propagators). Hence we can safely 
use the above propagator Eq. (j2J) with du > 3 (> 4) for scale (conformal) invariant 
analysis. As larger du values give smaller unparticle contributions, we would demand 
only scale invariance which allows smaller values of du- In addition, as the SM energy 
momentum tensor is a conserved quantity, it does not require any operator mixing under 
SM renormalisation. 

For our study, we closely followed Georgi's approach [7] to fix the coefficient of the 
effective interaction term at the weak scale for our phenomenology. Accordingly, the 
effective interaction at the weak scale is proportional to 

where Mu is the mass of heavy mediators in the hidden sector, Au is the scale at which 
dimensional transmutation occurs and Cu is the Wilson coefficient, dsz and du are the 
scaling dimensions of Banks-Zaks (BZ) [17] and the unparticle operators respectively. In 
the case of tensor unparticle operator coupled to the energy momentum tensor of the SM 
fields (having scaling dimension 4), the combination CuA.ff z /Mfc becomes dimensionless, 



which we call At in Eq. [TJ This can be phenomenologically interesting only if we assume 
the dimensional transmutation scale closer to the electroweak scale and hence we choose 
A u to be of the order of TeV. 

3 Diphotons at hadron colliders 

The hadronic cross section, da Pl,P2 , in the QCD improved parton model can be expressed 
in terms of mass factorized parton level subprocess cross sections, A ab convoluted with 
PDFs f a /p and f b/P : 

da Pl ' P2 (S) = ^2 J dXa J dx b f a/Pl (x a ,iJ,p) f b /p 2 (x b ,ijp) A ab (x a ,x b ,S,/j 2 F ) , (4) 

a,b=q,q,g 

where x a , x b are the momentum fractions of the incoming partons a and b respectively, 
S is the invariant mass square of the hadronic system and \x P is the factorization scale. 
Various kinematic constraints on the final states can be applied by introducing delta 
functions of these constraints. At LO in QCD we have contributions of order a 2 , a\ 2 , A*. 
Here a is the fine structure constant. At this order, the diphoton signal in the SM comes 
from the process q + q — > 7 + 7 and in the unparticle physics we have q + q — > 7 + 7 
and g + g — > 7 + 7. The presence of the gluon initiated process at LO in QCD is 
due to the fact that the unparticle fields couple to quarks and gluons with the same 
strength. At NLO, contributions are of the order a s a 2 , a s a\ 2 , and a s \f. In the SM, 
q + q — >7 + 7 + one loop, q + q — > 7 + 7 + g and q(q) + # — > 7 + 7 + q(q) contribute. 
In the unparticle case, in addition to the above processes, 5 , + g^7 + 7 + one loop and 
g + g — > 7 + 7 + fi', also contribute. In addition, at the NLO SM gg —>■ 77 box diagram 
contributes via its interference with leading order gg — > 77 diagram with an unparticle 
propagator. 

Before we compute the observable with the desired final state kinematical constraints, 
we compute the reduced cross sections (A ab in Eq. (j4])) that are free of both IR and 
ultraviolet singularities. The UV divergences resulting from the processes involving loops 
are regulated using dimensional regularization and renormalized in the MS scheme. 
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To study various kinematical distributions of these diphoton events with the ex- 
perimental cuts, a fully analytical computation at the NLO level is tedious. Hence 
we have opted for a semi-analytical approach, namely, a two cutoff phase-space slicing 
method [18]. Here, the phase space integrals are split into regions that are sensitive to 
soft and collinear singularities and those that are free of them. Such a slicing of the 
phase space is done using two very small parameters denoted by 5 S and 5 C which define 
the boundaries of soft and collinear regions respectively. Processes involving real gluons 
or quark/antiquark emissions can be decomposed as 

d<T l = daXsoft^s) + daZZtiS., S e ) + da^} m (5 s , S e ) . (5) 

The subscripts col and fin correspond to collinear and finite parts respectively. In the 
above equation, da^ j t (5 s ) is the cross section wherein the phase space integrals of the 
out going gluons are constrained such that the energy of the gluon in the CM frame of the 
incoming partons, is in the range < E g < 5 sv /s^/2. s a6 is the invariant mass squared 
of the incoming partons a and b. Similarly, da 1 ^ ol (S s ,S c ) is the cross section where 
the collinear regions of the final state partons are integrated out keeping the invariant 
mass squared of the pair of collinear partons to be less than 5 c s a b. The remaining cross 
section where the soft and collinear regions are removed is denoted by da^j in (S s ,S c ) 
and can be evaluated in 4-dimensions itself. The integrals involving soft and collinear 
regions of the phase space are regulated using dimensional regularization and can be 
computed analytically. The next step involves the computation of the virtual gluon 
corrections denoted by a\ to the born processes to NLO in QCD. The soft and collinear 
singularities present in the loop integrals of these virtual corrections are regulated using 
dimensional regularization. The soft and collinear singularities coming from the phase 
space as well as from the loop integrals manifest themselves as poles in e. One can easily 
show that the following combination 

da s ab +v (5 s ) = da:^ oft (5 s ) + aV, (6) 

is free of e poles coming from the soft gluons. The collinear singularities from the a\ h 
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and da^: ol (5 s , S c ) will go away if the mass factorization counter terms denoted by da F b , 
defined in the MS scheme, are included in the cross sections. The combination 

da s ab +v+c+F (5 s , 5 C ) = da s ab +v (5 s ) + do^tf., 5 C ) + da F ab , (7) 

is free of all IR singularities in QCD. Finally we end up with IR finite cross section given 
by 

^ = da s ab +v+c+F (5 s , S e ) + da2% n (6 s , S e ) . (8) 

Even though da^ b +v+c+F and daH a j in are separately dependent on 5 S and 5 C , their sum 
is expected to be independent of these parameters. Compact and simple analytical 
expressions for soft and collinear parts of the cross sections are reserved for a detailed 
publication [19]. 

The reduced cross section Eq. (jSJ) is IR finite. The additional final state QED singu- 
larities coming from the SM subprocess q(q) g — > 77 q(q) are suppressed by the smooth 
cone isolation criteria advocated by Frixione [10]. The isolation criteria on each of the 
photons is imposed by using the longitudinal boost invariant dimensionless parameter 
Rij = a/ (jji — rj-f) 2 + ((pi — 7 ) 2 where r\i and ^ 7 are pseudo-rapidities of the outgoing 
hadron i and the photon respectively. Similarly, <pi and 7 are their azimuthal angles 
with respect to the beam direction. Within a circle of radius Rq centered at each photon 
we reject any event unless the following condition is fulfilled: 

E r,i 0(R - Ri-y) < H(R) , for all R < R , (9) 

i 

where Eti is the transverse energy of the hadron i. The function 7i(R) must vanish as 
R goes to zero in order to get IR finite observable. One such choice put forth by Frixione 
is 

where E^° is a fixed energy and n > 1. 
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Since the isolation criteria on diphoton events does not allow any hard partonic 
(hadronic) activity closer (collinear) to either of the photons, the fragmentation contri- 
butions which are collinear in nature are suppressed. The isolation criteria also ensures 
the cancellation of the soft singularities coming from the real and virtual gluons as it 
does not restrict the phase space of the soft gluons. By imposing finite transverse mo- 
mentum cuts on the isolated photons, the smooth cone isolation criteria does not reject 
any hard parton which could be collinear to the initial states — ensuring mass factoriza- 
tion. Hence the isolation criteria on the diphotons guarantees the use of the phase space 
slicing method to compute an IR safe observable. In addition, we can safely employ the 
experimental cuts on the rapidity |y 7 | of the photons. 

4 Results 

In this section we present our results for = 14 TeV at the LHC. The virtual and 
real corrections have been evaluated in the limit of vanishing quark masses. We have 
evaluated LO and NLO cross sections with CTEQ6L and CTEQ6M parton density sets 
[20] with the corresponding strong coupling constant a s (M z ) = 0.118. The fine structure 
constant is a = 1/128. A single scale Q, the invariant mass of the photon pair, is used for 
the renormalization and factorization scale. The scale A u at which scale invariance sets 
in the BZ sector is chosen to be 1 TeV. Scale invariance restricts the scaling dimension of 
tensor operators to du > 3, in our calculation we have chosen du = 3.01. The coupling 
\ t is taken to be of order one. 

We have imposed the kinematical cuts on the photons as used by the ATLAS detector 
[21]: > 40 (25) GeV for harder (softer) photon, |?7 T | < 2.5 for each photon. The 
photons are restricted to have a separation of at least i? 77 = 0.4. In addition to these 
constraints, we impose smooth cone isolation criterion on photons given in Eq. (8). 
R = 0.4 is taken for the cone radius. Unless otherwise specified we use n = 2 and 
E l T so = 15 GeV which appear in H(R) 

To show that the results are independent of the choice of the slicing parameters 
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Figure 1: Plots showing stability of da/dQ for the SM (left panel) and the SM+ UP 
(right panel) against 5 S variation with the choice of S c = 10~ 5 . 

5 S and 8 C , we have plotted da^~ v+c+F (5 s ,5 c ) and da r a l a j in {S s , 5 C ) pieces of da/dQ as 
functions of 5 S with 5 C kept fixed at a very small value 10~ 5 in Fig. [TJ We see that 
the sum is fairly stable under the variation. The variation is around 2 % for SM and 
is less then 5 % for the signal (SM+UP). For rest of our numerical study, we have 
chosen 5 S = 1CT 3 and 5 C = 1CT 5 . As SM NLO results exist in literature [4-6], to further 
check our code we have compared our results with [5] using their isolation criterion 
[H(R) = Pr(l) e ([1 — cos(i?)]/[l — cos(-R )]) n ] an d their choice of ftp, /i^ and PDFs. We 
have found good agreement with [5]. 

We present various subprocess contributions to NLO in QCD in the invariant mass 
distribution for the range 400 < Q < 900 GeV in Fig. [2J In the SM both qq and qg 
subprocess contributions are positive with qq contribution being dominant over that of 
qg for the range of Q considered. However for smaller values of Q(< 150 GeV), qg 
contribution will be dominant due to the large gluon flux. In the unparticle sector the qq 
and gg subprocess contributions via the pure unparticle exchange (direct) are positive 
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Figure 2: Subprocess contributions in the SM and unparticle model in the invariant 
mass distribution at NLO for 400 < Q < 900 GeV. 



while the qg contribution is negative. At the interference level the gg interference with the 
SM box has a positive contribution and is almost constant for the range of Q considered. 
However, the interference of both qq and qg subprocesses with the SM have negative 
contributions and are larger in magnitude compared to the direct ones. As a result, the 
total W-particle contributions (direct+interference) can have a pattern that will change 
the sign at some high value of Q. Consequently, below this Q the signal can be lower 
than the SM background. The direct ^/-particle exchange contributions can become 
significant in the large Q region because the cross sections go as powers of Q/Ay, thus 
leading to the visibility of the W-particles only in that region of Q. It is worth noting 
that only the gg initiated subprocess has the dominant contribution over the rest in this 
region. 

In Fig. 3, we present our results for da/ dQ (left panel) for Q between 400 and 900 GeV 
and da/dY (right panel) for \Y\ < 2.0, where Y is the rapidity of the diphoton system. 
The short dashed (LO) and the long dashed lines (NLO) correspond to SM distributions. 
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Figure 3: Plots showing invariant mass (left panel) and rapidity (right panel) distribu- 
tions of the diphoton system with du = 3.01, A u = 1 TeV and X t = 0.9. For rapidity 
distribution Q is integrated in the range 600 GeV < Q < 0.9A W . 



The dotted (LO) and the solid (NLO) lines correspond to the signal (SM+UP) of the 
unparticle physics. As expected we find that the unparticle contribution to the invariant 
mass distribution grows with Q and it dominates over the SM contributions above Q = 
600 GeV. The precise value where this happens will depend very much on the choice of 
Au and other unparticle parameters. Since unparticle effects can be seen in the larger 
values of Q, rapidity distributions are computed by integrating Q between 600 and 900 
GeV. Near the central value of Y, we find large enhancement of the cross section from 
the SM results if we include unparticle contributions. 

At next-to-leading order the total transverse energy of the hadrons is due to a single 
parton around the photons and does not correspond to the actual hadronic energy in an 
experiment. Hence the E^° at the parton level is a crude estimate of that of the jets 
at the detector level. To show the dependence on we present the invariant mass 
distribution for = 15 GeV and = 30 GeV, for n = 2. To study the dependence 
of our predictions on the choice of 7i{R) we have varied it by changing n from 1 to 2 
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Figure 4: Dependence on isolation parameters (left panel) and n (right panel) is 
shown in invariant mass distribution for the range 400 < Q < 900 GeV. 

and keeping fixed at a value of 15 GeV. These variations are shown in Fig. [4] and 
show a very small dependence for R = 0.4. 

Until now our analysis was restricted to the case where 3 < du < 4, which was 
essential for a tensor unparticle as a consequence of scale invariance. There are no known 
examples of unitary quantum field theory that are scale invariant but not conformal 
invariant. For conformal invariance, unitarity demands that dy > 4 for the tensor 
unparticles. In Fig. (left panel) the unparticle sector as a result of scale invariance 
and not conformal (du = 3.01) is contrasted to the case where the unparticle sector is 
conformal (d u = 4.001). This is to both LO and NLO in QCD and for A u = 1 TeV. 
In Fig. [5] (right panel) we have considered the invariant mass distribution for du > 4 to 
LO in QCD. For this plot we have considered Au = 2 TeV and have probed Q < 0.9A^. 
Closer to du = 4 there could still be sufficient unparticle contribution for the tensor 
operator. The turn-around with energy behaviour of the unparticle effects is a typical 
feature of any physics beyond the SM. One observes a similar behaviour in models with 
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Figure 5: Invariant mass distribution of the diphoton system (left panel) for dy > 4 
(conformal invariance) is contrasted with the 3 < du < 4 (scale invariance). To LO in 
QCD and for A w = 2 TeV (right panel), we have plotted the invariant mass distribution 
of the diphoton system for du > 4. 

large extra-dimensions. The origin of this turn around contribution comes from the 
terms proportional to (s/Au) du in the matrix elements involving unparticles. Finally 
the observable we have considered for the diphoton process does not distinguish other 
beyond standard model contact interactions. 

5 Conclusions 



In this article we have used one of the important channels namely the diphoton produc- 
tion in order to explore the scale invariant unparticle sector coupled to the SM. We have 
restricted ourselves to contributions coming from the tensor unparticles. Since the QCD 
plays an important role at the LHC, we have systematically included all the partonic 
subprocesses that enter at the NLO level in the strong coupling constant a s . We have 
employed smooth cone isolation criteria advocated by Frixione in order to reduce frag- 
mentation contributions and thereby removing final state QED singularities. We have 
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applied appropriate cuts on the photon kinematics that are used by the ATLAS Collab- 
oration. Since it is hard to perform a fully analytic computation with isolation criteria 
and experimental cuts on the final state photons, we have resorted to the two cut-off 
phase space slicing method. We have also shown that our results are insensitive to the 
slicing parameters that bear no physics but are introduced in the intermediate stages of 
the computation. Our SM results are in agreement with the ones in the literature. We 
have presented two important distributions, namely, the invariant mass and the rapidity 
distributions of the diphoton system with the above isolation criteria and the relevant 
experimental cuts. For Au — 1 TeV, we have presented various subprocess contributions 
to the signal in the invariant mass distribution and found that the gg contribution is 
dominant at high Q values where the unparticle effects are visible. Both in the SM 
and in the unparticle case, the NLO contributions resulting from a s enhance the cross 
sections, thanks to the large gluon flux at the LHC. Our NLO results to the diphoton 
production including the unparticle effects not only estimate the corrections from QCD 
to the leading contributions but also reduce the scale uncertainties coming from the fac- 
torization scale [19]. We have also considered the case where the unparticle sector is a 
result of scale invariance and not conformal (3 < d u < 4) and also the case where it is 
due to conformal invariance {du > 4) for the tensor unparticle operator. 
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